Optomechanically Actuated Hydrogel Platform for Cell Stimulation with Spatial and Temporal Resolution.

Cells exist in the body in mechanically dynamic environments, yet the vast majority of in vitro cell culture is conducted on static materials such as plastic dishes and gels. To address this limitation, we report an approach to transition widely used hydrogels into mechanically active substrates by doping optomechanical actuator (OMA) nanoparticles within the polymer matrix. OMAs are composed of gold nanorods surrounded by a thermoresponsive polymer shell that rapidly collapses upon near-infrared (NIR) illumination. As a proof of concept, we crosslinked OMAs into laminin-gelatin hydrogels, generating up to 5 µm deformations triggered by NIR pulsing. This response was tunable by NIR intensity and OMA density within the gel and is generalizable to other hydrogel materials. Hydrogel mechanical stimulation enhanced myogenesis in C2C12 myoblasts as evidenced by ERK signaling, myocyte fusion, and sarcomeric myosin expression. We also demonstrate rescued differentiation in a chronic inflammation model as a result of mechanical stimulation. This work establishes OMA-actuated biomaterials as a powerful tool for in vitro mechanical manipulation with broad applications in the field of mechanobiology.

Mechanosensitive traction force generation is regulated by the neutrophil activation state.

The generation of traction forces by neutrophils regulates many crucial effector functions responsible for host defense, such as attachment, spreading, migration, phagocytosis, and NETosis. The activation state of the cell is a strong determinant of the functional efficacy of the neutrophil; however, the effect of activation on traction force production has not yet been determined experimentally. Previously, the mapping of cellular-generated forces produced by human neutrophils via a Traction Force Microscopy (TFM) method has required a three-dimensional imaging modality to capture out-of-plane forces, such as confocal or multiphoton techniques. A method newly developed in our laboratories can capture out-of-plane forces using only a two-dimensional imaging modality. This novel technique-combined with a topology-based single particle tracking algorithm and finite element method calculations-can construct high spatial frequency three-dimensional traction fields, allowing for traction forces in-plane and out-of-plane to the substrate to now be differentially visualized and quantified with a standard epifluorescence microscope. Here we apply this technology to determine the effect of neutrophil activation on force generation. Sepsis is a systemic inflammatory response that causes dysregulated neutrophil activation in vivo. We found that neutrophils from septic patients produced greater total forces than neutrophils from healthy donors and that the majority of this dysregulation occurred in-plane to the substrate. Ex vivo activation of neutrophils from healthy donors showed differential consequences depending on activation stimuli with mechanosensitive force decreases observed in some cases. These findings demonstrate the feasibility of epifluorescence-based microscopy in mapping traction forces to ask biologically significant questions regarding neutrophil function.

Beyond symptomatic diagnosis of mild traumatic brain injury.

It is commonly assumed that there is no brain injury if there are no noticeable symptoms following a head impact. There is growing evidence that traumatic brain injuries can occur with no outward symptoms and that the damage from these injuries can accumulate over time resulting in disease and impairment later in life. It is time to rethink the role that symptoms play in traumatic brain injury and adopt a quantitative understanding of brain health at the cellular level to improve the way we diagnose, prevent, and ultimately heal brain injury.

It is commonly assumed that there is no brain injury if there are no noticeable symptoms following a head impact. There is growing evidence that traumatic brain injuries can occur with no outward symptoms and that the damage from these injuries can accumulate over time resulting in disease and impairment later in life. It is time to rethink the role that symptoms play in traumatic brain injury and adopt a quantitative understanding of brain health at the cellular level to improve the way we diagnose, prevent, and ultimately heal brain injury.

Monocytes use protrusive forces to generate migration paths in viscoelastic collagen-based extracellular matrices.

Circulating monocytes are recruited to the tumor microenvironment, where they can differentiate into macrophages that mediate tumor progression. To reach the tumor microenvironment, monocytes must first extravasate and migrate through the type-1 collagen rich stromal matrix. The viscoelastic stromal matrix around tumors not only stiffens relative to normal stromal matrix, but often exhibits enhanced viscous characteristics, as indicated by a higher loss tangent or faster stress relaxation rate. Here, we studied how changes in matrix stiffness and viscoelasticity, impact the three-dimensional migration of monocytes through stromal-like matrices. Interpenetrating networks of type-1 collagen and alginate, which enable independent tunability of stiffness and stress relaxation over physiologically relevant ranges, were used as confining matrices for three-dimensional culture of monocytes. Increased stiffness and faster stress relaxation independently enhanced the 3D migration of monocytes. Migrating monocytes have an ellipsoidal or rounded wedge-like morphology, reminiscent of amoeboid migration, with accumulation of actin at the trailing edge. Matrix adhesions and Rho-mediated contractility were dispensable for monocyte migration in 3D, but migration did require actin polymerization and myosin contractility. Mechanistic studies indicate that actin polymerization at the leading edge generates protrusive forces that open a path for the monocytes to migrate through in the confining viscoelastic matrices. Taken together, our findings implicate matrix stiffness and stress relaxation as key mediators of monocyte migration and reveal how monocytes use pushing forces at the leading edge mediated by actin polymerization to generate migration paths in confining viscoelastic matrices. Significance Statement: Cell migration is essential for numerous biological processes in health and disease, including for immune cell trafficking. Monocyte immune cells migrate through extracellular matrix to the tumor microenvironment where they can play a role in regulating cancer progression. Increased extracellular matrix (ECM) stiffness and viscoelasticity have been implicated in cancer progression, but the impact of these changes in the ECM on monocyte migration remains unknown. Here, we find that increased ECM stiffness and viscoelasticity promote monocyte migration. Interestingly, we reveal a previously undescribed adhesion-independent mode of migration whereby monocytes generate a path to migrate through pushing forces at the leading edge. These findings help elucidate how changes in the tumor microenvironment impact monocyte trafficking and thereby disease progression.

Modeling high strain-rate microcavitation in soft materials: the role of material behavior in bubble dynamics.

Inertial Microcavitation Rheometry (IMR) is a promising tool for characterizing the mechanical behavior of soft materials at high strain rates. In IMR, an isolated, spherical microbubble is generated inside a soft material, using either a spatially-focused pulsed laser or focused ultrasound, to probe the mechanical behavior of the soft material at high strain rates (>103 s-1). Then, a theoretical modeling framework for inertial microcavitation, incorporating all the dominant physics, is used to extract information regarding the mechanical behavior of the soft material by fitting model predictions to the experimentally measured bubble dynamics. To model the cavitation dynamics, approaches based on extensions of the Rayleigh-Plesset equation are commonly used; however, these approaches cannot consider bubble dynamics that involves appreciable compressible behavior and place a limit on the nonlinear viscoelastic constitutive models that may be employed to describe the soft material. To circumvent these limitations, in this work, we develop a finite-element-based numerical simulation capability for inertial microcavitation of spherical bubbles that enables appreciable compressibility to be accounted for and more complex viscoelastic constitutive laws to be incorporated. We first apply the numerical simulation capability to understanding the role of material compressibility during violent spherical bubble collapse, and based on finite-element simulations, we propose a Mach-number-based threshold of 0.08, beyond which bubble collapse is violent, and the bubble dynamics involves compressibility not accounted for in Rayleigh-Plesset-based approaches. Second, we consider more complex viscoelastic constitutive models for the surrounding material, including nonlinear elastic and power-law viscous behavior, and apply IMR by fitting computational results to experimental data from inertial microcavitation of polyacrylamide (PA) gels in order to determine material parameters for PA gels at high strain rates.

Head Impact Modeling to Support a Rotational Combat Helmet Drop Test.

INTRODUCTION: The Advanced Combat Helmet (ACH) military specification (mil-spec) provides blunt impact acceleration criteria that must be met before use by the U.S. warfighter. The specification, which requires a helmeted magnesium Department of Transportation (DOT) headform to be dropped onto a steel hemispherical target, results in a translational headform impact response. Relative to translations, rotations of the head generate higher brain tissue strains. Excessive strain has been implicated as a mechanical stimulus leading to traumatic brain injury (TBI). We hypothesized that the linear constrained drop test method of the ACH specification underreports the potential for TBI. MATERIALS AND METHODS: To establish a baseline of translational acceleration time histories, we conducted linear constrained drop tests based on the ACH specification and then performed simulations of the same to verify agreement between experiment and simulation. We then produced a high-fidelity human head digital twin and verified that biological tissue responses matched experimental results. Next, we altered the ACH experimental configuration to use a helmeted Hybrid III headform, a freefall cradle, and an inclined anvil target. This new, modified configuration allowed both a translational and a rotational headform response. We applied this experimental rotation response to the skull of our human digital twin and compared brain deformation relative to the translational baseline. RESULTS: The modified configuration produced brain strains that were 4.3 times the brain strains from the linear constrained configuration. CONCLUSIONS: We provide a scientific basis to motivate revision of the ACH mil-spec to include a rotational component, which would enhance the test's relevance to TBI arising from severe head impacts.

Adhesive peptide and polymer density modulate 3D cell traction forces within synthetic hydrogels.

Cell-extracellular matrix forces provide pivotal signals regulating diverse physiological and pathological processes. Although mechanobiology has been widely studied in two-dimensional configurations, limited research has been conducted in three-dimensional (3D) systems due to the complex nature of mechanics and cellular behaviors. In this study, we established a platform integrating a well-defined synthetic hydrogel system (PEG-4MAL) with 3D traction force microscopy (TFM) methodologies to evaluate deformation and force responses within synthetic microenvironments, providing insights that are not tractable using biological matrices because of the interdependence of biochemical and biophysical properties and complex mechanics. We dissected the contributions of adhesive peptide density and polymer density, which determines hydrogel stiffness, to 3D force generation for fibroblasts. A critical threshold of adhesive peptide density at a constant matrix elasticity is required for cells to generate 3D forces. Furthermore, matrix displacements and strains decreased with matrix stiffness whereas stresses, and tractions increased with matrix stiffness until reaching constant values at higher stiffness values. Finally, Rho-kinase-dependent contractility and vinculin expression are required to generate significant 3D forces in both collagen and synthetic hydrogels. This research establishes a tunable platform for the study of mechanobiology and provides new insights into how cells sense and transmit forces in 3D.

Predicting complex nonspherical instability shapes of inertial cavitation bubbles in viscoelastic soft matter.

Inertial cavitation in soft matter is an important phenomenon featured in a wide array of biological and engineering processes. Recent advances in experimental, theoretical, and numerical techniques have provided access to a world full of nonlinear physics, yet most of our quantitative understanding to date has been centered on a spherically symmetric description of the cavitation process in water. However, cavitation bubble growth and collapse rarely occur in a perfectly symmetrical fashion, particularly in soft materials. Predicting the onset of dynamically arising, nonspherical instabilities in soft matter has remained a significant, unresolved challenge, in part due to the additional constitutive complexities introduced by the surrounding nonlinear viscoelastic solid. Here, we provide a new theoretical framework capable of accurately predicting the onset of nonspherical instability shapes of a bubble in a soft material by explicitly accounting for all pertinent nonlinear interactions between the cavitation bubble and the solid surroundings. Comparison with high-resolution experimental images from laser-induced cavitation events in a polyacrylamide hydrogel show excellent agreement. Interestingly, and consistent with experimental findings, our model predicts the emergence of various dynamic instability shapes for circumferential bubble stretch ratios greater than 1, in contrast to most quasistatic investigations. Our new theoretical framework not only provides unprecedented insight into the cavitation dynamics in a soft, nonlinear solid, but also provides a quantitative means of interpreting bubble dynamics relevant to a wide array of engineering and medical applications as well as natural phenomena.

Characterizing viscoelastic materials via ensemble-based data assimilation of bubble collapse observations.

Viscoelastic material properties at high strain rates are needed to model many biological and medical systems. Bubble cavitation can induce such strain rates, and the resulting bubble dynamics are sensitive to the material properties. Thus, in principle, these properties can be inferred via measurements of the bubble dynamics. Estrada et al. (2018) demonstrated such bubble-dynamic high-strain-rate rheometry by using least-squares shooting to minimize the difference between simulated and experimental bubble radius histories. We generalize their technique to account for additional uncertainties in the model, initial conditions, and material properties needed to uniquely simulate the bubble dynamics. Ensemble-based data assimilation minimizes the computational expense associated with the bubble cavitation model, providing a more efficient and scalable numerical framework for bubble-collapse rheometry. We test an ensemble Kalman filter (EnKF), an iterative ensemble Kalman smoother (IEnKS), and a hybrid ensemble-based 4D-Var method (En4D-Var) on synthetic data, assessing their estimations of the viscosity and shear modulus of a Kelvin-Voigt material. Results show that En4D-Var and IEnKS provide better moduli estimates than EnKF. Applying these methods to the experimental data of Estrada et al. (2018) yields similar material property estimates to those they obtained, but provides additional information about uncertainties. In particular, the En4D-Var yields lower viscosity estimates for some experiments, and the dynamic estimators reveal a potential mechanism that is unaccounted for in the model, whereby the apparent viscosity is reduced in some cases due to inelastic behavior, possibly in the form of material damage occurring at bubble collapse.

Breast tumor stiffness instructs bone metastasis via maintenance of mechanical conditioning.

While the immediate and transitory response of breast cancer cells to pathological stiffness in their native microenvironment has been well explored, it remains unclear how stiffness-induced phenotypes are maintained over time after cancer cell dissemination in vivo. Here, we show that fibrotic-like matrix stiffness promotes distinct metastatic phenotypes in cancer cells, which are preserved after transition to softer microenvironments, such as bone marrow. Using differential gene expression analysis of stiffness-responsive breast cancer cells, we establish a multigenic score of mechanical conditioning (MeCo) and find that it is associated with bone metastasis in patients with breast cancer. The maintenance of mechanical conditioning is regulated by RUNX2, an osteogenic transcription factor, established driver of bone metastasis, and mitotic bookmarker that preserves chromatin accessibility at target gene loci. Using genetic and functional approaches, we demonstrate that mechanical conditioning maintenance can be simulated, repressed, or extended, with corresponding changes in bone metastatic potential.

Acoustic cavitation rheometry.

Characterization of soft materials is challenging due to their high compliance and the strain-rate dependence of their mechanical properties. The inertial microcavitation-based high strain-rate rheometry (IMR) method [Estrada et al., J. Mech. Phys. Solids, 2018, 112, 291-317] combines laser-induced cavitation measurements with a model for the bubble dynamics to measure local properties of polyacrylamide hydrogel under high strain-rates from 103 to 108 s-1. While promising, laser-induced cavitation involves plasma formation and optical breakdown during nucleation, a process that could alter local material properties before measurements are obtained. In the present study, we extend the IMR method to another means to generate cavitation, namely high-amplitude focused ultrasound, and apply the resulting acoustic-cavitation-based IMR to characterize the mechanical properties of agarose hydrogels. Material properties including viscosity, elastic constants, and a stress-free bubble radius are inferred from bubble radius histories in 0.3% and 1% agarose gels. An ensemble-based data assimilation is used to further help interpret the obtained estimates. The resulting parameter distributions are consistent with available measurements of agarose gel properties and with expected trends related to gel concentration and high strain-rate loading. Our findings demonstrate the utility of applying IMR and data assimilation methods with single-bubble acoustic cavitation data for measurement of viscoelastic properties.

Epifluorescence-based three-dimensional traction force microscopy.

We introduce a novel method to compute three-dimensional (3D) displacements and both in-plane and out-of-plane tractions on nominally planar transparent materials using standard epifluorescence microscopy. Despite the importance of out-of-plane components to fully understanding cell behavior, epifluorescence images are generally not used for 3D traction force microscopy (TFM) experiments due to limitations in spatial resolution and measuring out-of-plane motion. To extend an epifluorescence-based technique to 3D, we employ a topology-based single particle tracking algorithm to reconstruct high spatial-frequency 3D motion fields from densely seeded single-particle layer images. Using an open-source finite element (FE) based solver, we then compute the 3D full-field stress and strain and surface traction fields. We demonstrate this technique by measuring tractions generated by both single human neutrophils and multicellular monolayers of Madin-Darby canine kidney cells, highlighting its acuity in reconstructing both individual and collective cellular tractions. In summary, this represents a new, easily accessible method for calculating fully three-dimensional displacement and 3D surface tractions at high spatial frequency from epifluorescence images. We released and support the complete technique as a free and open-source code package.

Flagellar kinematics reveals the role of environment in shaping sperm motility.

Swimming spermatozoa from diverse organisms often have very similar morphologies, yet different motilities as a result of differences in the flagellar waveforms used for propulsion. The origin of these differences has remained largely unknown. Using high-speed video microscopy and mathematical analysis of flagellar shape dynamics, we quantitatively compare sperm flagellar waveforms from marine invertebrates to humans by means of a novel phylokinematic tree. This new approach revealed that genetically dissimilar sperm can exhibit strikingly similar flagellar waveforms and identifies two dominant flagellar waveforms among the deuterostomes studied here, corresponding to internal and external fertilizers. The phylokinematic tree shows marked discordance from the phylogenetic tree, indicating that physical properties of the fluid environment, more than genetic relatedness, act as an important selective pressure in shaping the evolution of sperm motility. More broadly, this work provides a physical axis to complement morphological and genetic studies to understand evolutionary relationships.

Generating Cell Type-Specific Protein Signatures from Non-symptomatic and Diseased Tissues.

Here we demonstrate a technique to generate proteomic signatures of specific cell types within heterogeneous populations. While our method is broadly applicable across biological systems, we have limited the current work to study neural cell types isolated from human, post-mortem Alzheimer's disease (AD) and aged-matched non-symptomatic (NS) brains. Motivating the need for this tool, we conducted an initial meta-analysis of current, human AD proteomics studies. While the results broadly corroborated major neurodegenerative disease hypotheses, cell type-specific predictions were limited. By adapting our Formaldehyde-fixed Intracellular Target-Sorted Antigen Retrieval (FITSAR) method for proteomics and applying this technique to characterize AD and NS brains, we generated enriched neuron and astrocyte proteomic profiles for a sample set of donors (available at www.fitsarpro.appspot.com ). Results showed the feasibility for using FITSAR to evaluate cell-type specific hypotheses. Our overall methodological approach provides an accessible platform to determine protein presence in specific cell types and emphasizes the need for protein-compatible techniques to resolve systems complicated by cellular heterogeneity.

Application of mild hypothermia successfully mitigates neural injury in a 3D in-vitro model of traumatic brain injury.

Therapeutic hypothermia (TH) is an attractive target for mild traumatic brain injury (mTBI) treatment, yet significant gaps in our mechanistic understanding of TH, especially at the cellular level, remain and need to be addressed for significant forward progress to be made. Using a recently-established 3D in-vitro neural hydrogel model for mTBI we investigated the efficacy of TH after compressive impact injury and established critical treatment parameters including target cooling temperature, and time windows for application and maintenance of TH. Across four temperatures evaluated (31.5, 33, 35, and 37°C), 33°C was found to be most neuroprotective after 24 and 48 hours post-injury. Assessment of TH administration onset time and duration showed that TH should be administered within 4 hours post-injury and be maintained for at least 6 hours for achieving maximum viability. Cellular imaging showed TH reduced the percentage of cells positive for caspases 3/7 and increased the expression of calpastatin, an endogenous neuroprotectant. These findings provide significant new insight into the biological parameter space that renders TH effective in mitigating the deleterious effects of cellular mTBI and provides a quantitative foundation for the future development of animal and preclinical treatment protocols.

Mechanophenotyping of 3D multicellular clusters using displacement arrays of rendered tractions.

Epithelial tissues mechanically deform the surrounding extracellular matrix during embryonic development, wound repair, and tumor invasion. Ex vivo measurements of such multicellular tractions within three-dimensional (3D) biomaterials could elucidate collective dissemination during disease progression and enable preclinical testing of targeted antimigration therapies. However, past 3D traction measurements have been low throughput due to the challenges of imaging and analyzing information-rich 3D material deformations. Here, we demonstrate a method to profile multicellular clusters in a 96-well-plate format based on spatially heterogeneous contractile, protrusive, and circumferential tractions. As a case study, we profile multicellular clusters across varying states of the epithelial-mesenchymal transition, revealing a successive loss of protrusive and circumferential tractions, as well as the formation of localized contractile tractions with elongated cluster morphologies. These cluster phenotypes were biochemically perturbed by using drugs, biasing toward traction signatures of different epithelial or mesenchymal states. This higher-throughput analysis is promising to systematically interrogate and perturb aberrant mechanobiology, which could be utilized with human-patient samples to guide personalized therapies.

Context-Dependent Role of Vinculin in Neutrophil Adhesion, Motility and Trafficking.

Neutrophils are innate immune effector cells that traffic from the circulation to extravascular sites of inflammation. ß2 integrins are important mediators of the processes involved in neutrophil recruitment. Although neutrophils express the cytoskeletal protein vinculin, they do not form mature focal adhesions. Here, we characterize the role of vinculin in ß2 integrin-dependent neutrophil adhesion, migration, mechanosensing, and recruitment. We observe that knockout of vinculin attenuates, but does not completely abrogate, neutrophil adhesion, spreading, and crawling under static conditions. However, we also found that vinculin deficiency does not affect these behaviors in the presence of forces from fluid flow. In addition, we identify a role for vinculin in mechanosensing, as vinculin-deficient neutrophils exhibit attenuated spreading on stiff, but not soft, substrates. Consistent with these findings, we observe that in vivo neutrophil recruitment into the inflamed peritoneum of mice remains intact in the absence of vinculin. Together, these data suggest that while vinculin regulates some aspects of neutrophil adhesion and spreading, it may be dispensable for ß2 integrin-dependent neutrophil recruitment in vivo.

Modeling tissue-selective cavitation damage.

The destructive growth and collapse of cavitation bubbles are used for therapeutic purposes in focused ultrasound procedures and can contribute to tissue damage in traumatic injuries. Histotripsy is a focused ultrasound procedure that relies on controlled cavitation to homogenize soft tissue. Experimental studies of histotripsy cavitation have shown that the extent of ablation in different tissues depends on tissue mechanical properties and waveform parameters. Variable tissue susceptibility to the large stresses, strains, and strain rates developed by cavitation bubbles has been suggested as a basis for localized liver tumor treatments that spare large vessels and bile ducts. However, field quantities developed within microns of cavitation bubbles are too localized and transient to measure in experiments. Previous numerical studies have attempted to circumvent this challenge but made limited use of realistic tissue property data. In this study, numerical simulations are used to calculate stress, strain, and strain rate fields produced by bubble oscillation under histotripsy forcing in a variety of tissues with literature-sourced viscoelastic and acoustic properties. Strain field calculations are then used to predict a theoretical damage radius using tissue ultimate strain data. Simulation results support the hypothesis that differential tissue responses could be used to design tissue-selective treatments. Results agree with studies correlating tissue ultimate fractional strain with resistance to histotripsy ablation and are also consistent with experiments demonstrating smaller lesion size under exposure to higher frequency waveforms. Methods presented in this study provide an approach for modeling tissue-selective cavitation damage in general.

Comparative study of the dynamics of laser and acoustically generated bubbles in viscoelastic media.

Experimental observations of the growth and collapse of acoustically and laser-nucleated single bubbles in water and agarose gels of varying stiffness are presented. The maximum radii of generated bubbles decreased as the stiffness of the media increased for both nucleation modalities, but the maximum radii of laser-nucleated bubbles decreased more rapidly than acoustically nucleated bubbles as the gel stiffness increased. For water and low stiffness gels, the collapse times were well predicted by a Rayleigh cavity, but bubbles collapsed faster than predicted in the higher stiffness gels. The growth and collapse phases occurred symmetrically (in time) about the maximum radius in water but not in gels, where the duration of the growth phase decreased more than the collapse phase as gel stiffness increased. Numerical simulations of the bubble dynamics in viscoelastic media showed varying degrees of success in accurately predicting the observations.